Preprint typeset in JHEP style - HYPER VERSION hep-th/0406002 



(N 



O 



A dialogue of multipoles: matched asymptotic 
expansion for caged black holes 



o 
o 

(N 
^ ' Dan Gorbonos and Barak Kol 



Racah Institute of Physics, Hebrew University 
Jerusalem 91904, Israel 



CO ■ E-mail: gdanOphys . huj i . ac . il , barak_kol(Sphys . huj i . ac . il 

> 

(N 

O . 

^) , Abstract: No analytic solution is known to date for a black hole in a compact dinien- 

^^ \ sion. We develop an analytic perturbation theory where the small parameter is the size 

Tjj- I of the black hole relative to the size of the compact dimension. We set up a general pro- 

cedure for an arbitrary order in the perturbation series based on an asymptotic matched 
expansion between two coordinate patches: the near horizon zone and the asymptotic zone. 
ji» ■ The procedure is ordinary perturbation expansion in each zone, where additionally some 

D , boundary data comes from the other zone, and so the procedure alternates between the 

zones. It can be viewed as a dialogue of multipoles where the black hole changes its shape 
(mass multipoles) in response to the field (multipoles) created by its periodic "mirrors", 
K> ' and that in turn changes its field and so on. We present the leading correction to the 

C^ ' full metric including the first correction to the area-temperature relation, the leading term 

for black hole eccentricity and the "Archimedes effect". The next order corrections will 
appear in a sequel. On the way we determine independently the static perturbations of 
the Schwarzschild black hole in dimension d > 5, where the system of equations can be 
reduced to "a master equation" — a single ordinary differential equation. The solutions 
are hypergeometric functions which in some cases reduce to polynomials. 
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1. Introduction and summary 

1.1 Background 

In the presence of extra compact dimensions, there exist several phases of black objects, 
namely massive solutions of General Relativity, depending on the relative size of the object 
and the relevant length scales in the compact dimensions. For concreteness, we consider a 
background with a single compact dimension — ]R°'~^'^ x S^, where d is the total space- 
time dimension and d > 5 (in order to avoid spacetimes with 2 or less extended spatial 
dimensions where the presence of a massive source is inconsistent with asymptotic flatness). 

In this background one expects at least two phases of black object solutions: when the 
size of the black object is small (compared to the size of the extra dimension) one expects 
the region near the object to closely resemble a d-dimensional black hole, while as one 
increases the mass one expects that at some point the black hole will no longer fit in the 
compact dimension and a black string, whose horizon winds around the compact dimension 
will be formed. Thus we distinguish between the black hole and the black string according 
to their horizon topology which is either spherical — S or cylindrical — S x S^, 
respectively. We shall sometimes refer to such a black hole in a compact dimension as a 
"caged black hole". 

More generally one could consider backgrounds R'^~'^~^'^ x X*^ where X"^ is any c- 
dimensional compact Ricci-flat manifold such as the c-dimensional torus, the K3 surface 
or a Calabi-Yau 3-fold. For a general X more black object phases will exist, but we expect 
that generically the essential phase transition physics between any two specific phases will 
be qualitatively similar to the X = S^ case. 

This system raises several deep questions in general relativity [Q] : during the real-time 
phase transition a naked singularity may be encountered which may require intervention 
from quantum gravity and an amendment to Cosmic Censorship; the transition would cer- 
tainly produce some sort of a "thunderbolt" [^, ^; the system exhibits a critical dimension 
for stability in at least two instances f^, Q; it is a prototype example for the failure of 
black hole uniqueness in higher dimensions d > 5 |5|, ^; a novel kind of topology change 
is expected to play a central role [Q]; and there is an ongoing debate regarding the correct 
phase diagram, especially whether a stable non-uniform^ string phase exists (for all d) Q. 

Considerable work, much of it recent, went into finding the various solutions in this 
background. A uniform black string is readily described analytically by an arbitrary mass 
Schwarzschild solution in d — 1 dimensions with an added spectating coordinate. Gregory 
and Laflamme (GL,1994) discovered that this solution develops a tachyon below a certain 
critical mass M. Horowitz and Maeda (2000) gave an argument for the existence of sta- 
ble non-uniform strings [0|. Gubser analytically perturbed the critical GL string to find 
approximate expressions for non- uniform solutions (which we interpret to be unstable) Q . 
De-Smet attempted to find analytic solutions by classifying 5d algebraically special metrics 
(the same method Kerr used successfully in 4d to obtain the rotating black hole) not find- 
ing any novel ones in this background |1Q]. Harmark and Obers found a clever ansatz that 



^not invariant under translation along the compact direction. 



reduces the number of unknown metric functions from 3 to 2, but could not completely 
solve the equations either [0]. Recently Sorkin generalized the analysis of Gubser from 
5d (and 6d which was done by Wiseman ||l2|, fixing some problems with the earlier anal- 
ysis) to an arbitrary dimension, discovering an interesting critical dimension: for d > 14 
the non-uniform branch emanating from the GL point changes its nature and essentially 
becomes stable Q. A critical dimension for a different point in the phase diagram was 
already predicted in [|[. 

The limited success of analytical methods created a demand for numerical solutions. 
The branch of non-uniform black string solutions emanating from the GL point was ob- 



tained numerically by Wiseman |^2[ (see also earlier work in [13| and a post-analysis 
in IjlJ, |l^) who managed to formulate axially-symmetric gravitostatics (namely, essen- 
tially 2d) in a "relaxation" form (a procedure familiar from electrostatics) while presenting 
the constraints through "Cauchy-Riemann — like" relations. Even though there is no 
definitive answer yet whether these solutions are indeed unstable, the author argued that 
irrespectively they cannot serve as an endpoint for the decay of the GL string. While 
all the solutions we mentioned above are static Choptuik et al. performed a demanding 
numerical time evolution for the decay of the black string, but had to stop before the end 
state was reached due to an essential limitation of the algorithm used (grid stretching) in 



the high curvature region which forms |16] 



More recently the focus shifted from black strings to black holes. While intuition leads 
us to expect that a small black hole should exist being indifferent to the existence of a much 
larger compact dimension, no analytic solution is available to date. In [jl^ indications for 
the nature of the phase transition were gained from an analysis of possible time-symmetric 
initial data. In [fq] the closely related problem of black holes in a braneworld was tackled 



numerically. In [19, |2^ it was shown that indeed there are order parameters such that the 
black hole and black string are at finite values, as was assumed in Q, and moreover pO| 
announced most of the quantitative results of the current paper. In |21, 22 1 numerical black 



hole solutions were presented for the first time in 5d and 6d respectively, giving strong 
evidence for their existence. Finally, p3] presented a "first order analytic approximation" 



of small black holes in the framework of the Harmark-Obers coordinates |11|, and |24| 
found an analytic approximation for a small black hole on a brane. 

1.2 Motivation and basic set-up 

In this paper we present the first analytic (though perturbative) procedure to obtain so- 
lutions for small black holes (BH's). Let us introduce some notation (see figure ||). We 
denote by z, r the "cylindrical" coordinates, where z is the coordinate along the compact 
dimension whose period we denote by L, and r is the radial coordinate in the extended 
M spatial dimensions. The problem is characterized by a single dimensionless parame- 
ter, for instance the dimensionless mass /i = G^M/L'^^^ where Gn is the d-dimensional 
Newton constant and M is the mass (measured at infinity), or fip = f3/L where (3 is the 
inverse temperature. In the vicinity of the black hole it is useful to introduce "spherical" 
coordinates p, x S'S well. We denote by /?o the Schwarzschild radius of the BH (in the small 
BH limit), and one has Pq~ = const Gtv M where const is a dimensionless constant. 




Figure 1: Illustration of the (r, z) "cylindrical" coordinates and the (p, x) "spherical" near horizon 
coordinates. The period of the compact dimension (in the z direction) is denoted by L. 

There is reason to expect good analytic control of small black holes even if we do not 
have a complete analytic solution since we have two good approximation in two different 
regions which overlap: for p <C 1/ the metric is expected to resemble closely a d-dimensional 
Schwarzschild-Tangherlini BH ||2^ , while for p ^ po the gravitational field is weak and the 
newtonian approximation holds. Hence // or more precisely po/L ex p^' is our small 
parameter for the perturbation. 

The motivations for this research are first to obtain a theoretical description of this 
simple system which is important on its own right, and second to gain understanding of the 
phase transition physics through combination with numerical work. The symbiosis with 
numerical work comes close to serve as a partial substitute of experiments (which are sorely 
absent in this field): the numerics are essential for understanding big black holes close to 
the phase transition where the perturbative expansion is expected to break down, and the 
analytic control serves to formulate the aims and methods of the numerics. Moreover, 
the two can be used to test and confirm one another, as was the case for this research 
and [£^, 21]. As it turns out the largest BH's obtained numerically show only a single 

22[, and that lends some hope 



multipole mode correction to their spherical horizon |21 



that the analytic expansion would retain some validity for large BH's as well. 

Basic set-up. The first decision to be made it to choose the coordinates and ansatz 
for the analysis. At first one would hope to use a single coordinate patch for the whole 



metric. However, |18] showed that in the popular conformal coordinates (where the metric 
in the {r, z) plane is in conformal form ds'^ = e^ ''''^•' (dr^ + dz^), see also below ( 5.19| )) 
the coordinate size of the horizon is a conformal invariant and hence the coordinate patch 
necessarily changes with p. A similar phenomenon happens in the Harmark-Obers coordi- 
nates which are a semi-infinite cylinder with a single marked point [(r, z) = (0, L/2)] where 
the coordinate transformation is singular and whose location changes with p. 

Therefore we choose to work with two coordinate patches (see figure pi) : the near zone 
p <^ L where the horizon (po) is fixed and the periodicity of z is invisible far away, and 
the asymptotic zone p 3> po where L is fixed and po is invisible. The metric in the two 
regions must be consistent over the overlap region po ^ p ^ L (which grows indefinitely 
as p ^ 0). 
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Figure 2: The division of the spacetime into two zones: the near zone p <^ L where po is fixed 
and the perturbative parameter is L~^, and the asymptotic zone p ^ po where L is fixed and 
the perturbative parameter is po- The two zones overlap over the overlap region which increases 
indefinitely in the small black hole limit. During the perturbation process the two zones are separate, 
and communicate only through the matching "dialogue". The near zone is defined by {{p,x) '■ P^ 
Po} while the asymptotic zone is defined by {(r, z) : r > 0, z ^ z + i}\(0, 0). 



Such a procedure is known in General Relativity as a ''matched asymptotic expansion" 
the metric is solved for in each asymptotic region and certain quantities are determined 



by matching the metrics over the overlap (for some recent examples see [26, 27, 28 1, in 
mathematical physics this idea goes back as far as Laplace who used it to find the shape 
of a drop of liquid on a surface — see [^9[ and references therein for a historical review). 
However, this is probably the first time such a procedure is used to find a static black hole 
solution. 

We start by defining the domain for each zone and the zeroth order solution. In the 
near zone, whose domain is {(/>, x) '■ P^ Po} we have a d-dimensional Schwarzschild black 
hole metric, with fixed po and the perturbation is in orders of L~^. In the asymptotic zone, 
on the other hand, the domain is {(r, z) : r > 0, z ~ z + L}\{0, 0), the zeroth order solution 
is simply the flat "cylinder" with the origin omitted, and the perturbation parameter is po- 

Our objective is to describe the perturbation process (to any order) and apply it. The 
paper is organized as follows. In section |2| we describe the equations in the asymptotic 
zone and especially the newtonian potential. In section ^ we determine the linearized 
corrections to the Schwarzschild black hole in the near zone. In section ^ we describe the 
general perturbation procedure for this system and in |5| we present quantitative results 
on the leading corrections to the zeroth order metric. In the appendices we review some 
information on Heun's equation and the Hyper geometric equation, review the definition of 
the surface gravity and give some details on vector harmonics in 5d (on S'^). We now turn 
to the summary. 



1.3 Summary 

The asymptotic zone. In section |2| we describe the asymptotic zone. The first correction 
to the zeroth order metrics described above is readily computed — it is the newtonian 



approximation in the asymptotic region (in standard harmonic gauge) , where the newtonian 



potential (2.6) is obtained through the method of mirror images (considering the infinite 
sequence of sources in the covering space at z = nL for any integer n). This term is 
proportional to Gn M cc Pq ~ and hence belongs to order d — 3 in the asymptotic zone. 
Actually, the whole post-newtonian procedure is relevant and we review it, though in this 
paper all we need is the lowest order (newtonian) approximation. 

Black hole perturbations. The next correction to consider is the leading correction to 
the Schwarzschild solution. It is the analogue of the newtonian approximation only here 
the unperturbed background is curved, and it describes the response of the geometry to 
the mirror sources far away. Despite the analogy with the newtonian approximation the 
implementation is involved and is described in section ^. 

In 4d this computation was carried out by Regge and Wheeler [30| who solved for all the 



linear perturbations, not only the static ones. Recently Ishibashi and Kodama succeeded 
to generalize their result to an arbitrary dimension, and to include a cosmological constant 
as well [^, ^, ^. Our treatment is independent and we compare the two approaches 
below after describing our own. 

The computation involves a few steps. We start by writing down the most general 
static perturbation to the metric. The spherical symmetry guarantees that at linear order 
perturbations in different representations of the rotation group will not mix, and we find by 
counting degrees of freedom that it suffices to consider "scalar harmonics" — representa- 
tions which are the symmetric product of the vector representation, or equivalently, metric 
perturbations which are determined by scalar functions on the sphere. It turns out that the 
spherical symmetry also suggests a natural gauge which we term "no derivatives gauge" 
and completely fixes the reparameterization invariance, leaving us with 3 undetermined 
metric functions (fields). 

Writing down the equations of motion and separating the angular variables we find 
that a Ricci fiatness condition in the angular directions yields an algebraic relation among 
the radial functions (3.27) which is similar to a trace condition and allows us to eliminate 



one of the fields. After substitution one can express one of the remaining fields in terms 
of the other and its first and second derivatives. Performing the second substitution we 
are left with a second order ordinary differential equation (ODE), rather than a third or 
fourth order one would initially expect. So finally one has a single second order ODE in 
the radial direction, for one metric function (and for each spherical harmonic mode) from 
which the whole metric may be recovered. This is the so called "master equation" which 
after a change of variables simplifies further to become ( 3.43| ). It would be nice to have a 



deeper understanding why these reductions were to be expected. 

The master equation belongs to the Heun class of Fuchsian equations, where Fuchsian 
means that the equation has only regular-singular points on the complex sphere which 
includes infinity, and Heun means that there are exactly 4 such points. Unlike the Hyper- 
geometric case of 3 regular singularities there is no general solution to the Heun equation, 
though several methods are available. In this case however, it turns out that the solutions 
can be written in terms of a hypergeometric function, and it would be nice to understand 



why that had to be the case. Interestingly, we observe that in some of the relevant cases 
these hyper geometric functions simplify further to polynomials, and in particular in 5d all 
relevant solutions are polynomials (solutions which are of even multipole number and are 
regular at the horizon). 

We started working out this problem before we were aware of the results of Kodama 
and Ishibashi |31, 32, El] and we continued independently even after learning about these 



papers in order to avoid the formalism of gauge invariant perturbation theory, and the 
various changes of variables which are employed there. We were able to do so and actually 
found a somewhat different master (Heun) equation. Yet the final reduction of our master 
Heun equation to a hypergeometric one was motivated by those papers. 

The matching procedure. One of the main results of this paper is the construc- 
tion of a perturbation method for the metric (in both patches) which may be carried 
in principle to an arbitrarily high order in the small parameter. The method is de- 
scribed in section ^. A crucial step is to identify a dimensionful expansion parameter 
on each patch: pQ in the asymptotic zone and L~^ in the near region. As in any per- 
turbative expansion, at each order one needs to solve a non-homogenous linear equa- 
tion — the linear equation being the same as the one which appears at first order and 
the non-homogeneous source term being constructed from lower order metric functions. 
The precise form of the source term depends on the higher order gauge choice which 
we do not specify, but will not change the method we describe. The solution to this 
equation is determined up to a solution of the homogeneous equation. This indeter- 
minacy for each zone on its own reflects the freedom of adding external field multi- 
poles — in the asymptotic zone they are situated at the origin while in the near zone 
they are at infinity. These external multipoles must be determined by matching with 
the other zone, a procedure which requires to identify (after matching the gauge) the 
leading terms in the metric on both zones. We call this process "a dialogue of multi- 
poles". 

A priori it is not obvious that the required terms from the other zone are already 
available at the right time, namely that the method is well-posed (that there are sufficient 
boundary conditions). Hence it is interesting to study at any given order in a specific 
zone which orders must be already available for matching from the other one, and thereby 
describe the pattern of the dialogue — the orders at which one should alternate between the 
zones. This pattern can be determined by a simple dimensional analysis of the multipole 
coefficients as we describe in the text, and indeed we find the system to be well-posed. 
Interestingly, the dialogue pattern which emanates is dimension dependent: 5d is special in 
that one scales a single order in the perturbation ladder on each zone and then alternates 
to the other zone; for d > 5 one needs to climb several steps before going to the other zone, 
and in the d ^ oo limit one gets infinitely many constants already from matching with the 
newtonian potential alone. 

Quantitative matching results. In section ^ we apply the general procedure to the 
leading order and match the newtonian potential at the asymptotic zone to get the leading 
correction to the Schwarzschild solution in the near zone. To this purpose it is essential to 



have available certain matching constants which can be read from the explicit solutions we 
derived for the linear perturbations of Schwarzschild. The next order correction is currently 
under study |34]. 

^From the metric which we obtain one may extract certain "measurables" : 

• The leading correction to the mass — temperature relation is given in (^.Sj). At this 
order the BH is still spherical but there is a correction to this relation since the small 
black hole does not asymptote in the near zone to flat space with zero potential, but 
rather there is a non-zero potential shift due to the images. 

• The leading (quadrupole) departure from a spherical horizon — measured by the 



"eccentricity" is given in ( p. 181) . The result of the deformation is to make the black 
hole longer along the z axis compared to the r axis as in figure |9| and can be understood 
from the shape of small (newtonian) equipotential lines around p = (see figure ^. 



The coefficient of the "inter-polar distance" is given in ( 5.3C1| ). By "inter-polar" 



distance we mean the proper distance from the "north pole" of the black hole around 



the compact circle and up to the "south pole" (see figure 11). Actually the black 
hole tends to "make room" for itself, in the sense that the inter-polar distance added 
to the black hole size in conformal coordinates is always larger than L, the size of 
the compact dimension. This can be re-stated as the observation that such black 
holes seem to always have a positive scalar charge as seen from infinity similar to the 
ordinary positive mass theorem (where the scalar is the one which arises from the 
dimensional reduction of the gzz metric component — the size of the extra dimension) . 
In 5d the effect is the strongest, where to leading order in the small parameter the 
inter-polar distance does not decrease at all. We term that "a black hole Archimedes 
effect" since the black hole repels or expands an amount of space equal to its size in 
5d (and less in higher dimensions). 



Most of these results were already announced in |20| and here we add the determina- 



tion of the inter-polar distance and the generalization of the eccentricity for d > 5. They 



were numerically confirmed in 5d |2l[ as well as in 6d [pq, SC] and other dimensions pq 



Recently a paper |23| has appeared deriving the leading order form of the metric within the 



framework of the Harmark-Obers coordinates |11[|, and as such overlaps with the results 



announced in |20| and proven here. The overlap includes the corrections to the temper- 



ature and area, while [23| obtains also the corrections to the mass and tension, and this 
paper derives the eccentricity and "Archimedes effect". Moreover, here we go beyond and 
demonstrate a method for an arbitrary number of successive approximations. 



Note added: In the 3rd version we performed some rewriting of subsection 3.3.2 where 
we made minor corrections and improved the clarity and we corrected a trivial factor of 2 
(15331) . 



m 



2. The asymptotic zone 

In this section we write the static Einstein equations in a form that will be convenient 
for iterative expansion in a small parameter around the flat Minkowsky spacetime. This 



type of expansion is known as "post-newtonian expansion" (see [£9[) and we follow here 
the usual conventions for this type of expansion. The leading order in the expansion is the 
newtonian approximation. We repeat here the calculation of the newtonian approximation 
for the caged black hole which appeared in many places (see for example Jill , p3|). The 
following orders in the expansion are called "post-newtonian" and their calculation will 
appear in |34]. 



For the post-newtonian expansion it is convenient to write the Ricci tensor in the 



following form |37] 



1 rflnP'^ 

where T^ q,^ and F^ are the Christoffel symbols of the first and the second kind, respec- 
tively, and in addition one defines 

r =lf ^M ^_ ^r" 

i M- - 2 V^^" dx- + ^"" dx^ ^'^"^'^^ dx^ 

Next one chooses the harmonic (or de Bonder ^) gauge by the requirement that 

r- = Ux^ = -^^{,/^ggP-) = 0. (2.2) 

^J—g oxP 

where we denote by g the determinant of the metric g^^. In this gauge, the last term in 
the expression of the Ricci tensor above vanishes. This choice of gauge is very convenient 
for expansion in the asymptotic zone. Finally one attempts to solve Einstein's equations. 
The first step in this iterative procedure is to look at the linearized equations valid for 
weakly gravitating regions, namely making the newtonian approximation. The metric is 
taken to be 

g^lv = rji^u + h/iu ■ (2.3) 

The harmonic gauge equation (^^) takes the more famous form (for example in the treat- 
ment of gravitational waves) 

Ku-^h^v'"'^ =0. 



One defines 



hf,u = h^,y--h'^r]f,^, (2.4) 



in terms of which the linearized field equations become 



^nV = G^, = 8^GdT^,, (2.5) 



where D is the flat space D'alambertian. 
^The first introduction of tliis gauge appeared in 



In our case, working on the covering space implies an infinite array of newtonian sources 
in the z direction and the only non-zero components of the energy momentum tensor is 

oo 

Too= Y. M6''-\xi,...,Xd-2)S{z-nL), 

n=— oo 

where (xi, . . . ,Xd-2) denote the extended spatial coordinates. The method of images can 
be used to solve the equation for hoo 

- _ $7V _ lOvrGrfM ^ 1 

^00 - -7- = jj^^TTw— 2^ — — : -7-1=^ 



d-2., - 1 



2 



where $Ar is the newtonian potential, conventionally normalized such that its flux through 
a surface enclosing a mass M is AirGiiM, pQ is given by |3^ 

.3 _ IGTTGdM 
"' - {d- 2)^^-2' ^ ' 

and 

d-l 

is the area of a unit S* . A more formal alternative to obtain the prefactor of the 
newtonian potential in equation ( |2.6D is through matching with the Schwarzschild metric 
in the near zone. 

We see that the first correction to the metric in the asymptotic region is of order p^ ~ . 
We will see later that we must choose the perturbation parameter in this region to be po 

rather than p^ ~ , namely^ 

00 

ff/.. = E^o"5&^ (2.8) 

n=0 

and hence we see that the leading correction comes at order d — 3, namely 

5it"'^=V- (2-9) 

Transforming back to /i^;/ using the inverse of ( |2.4| ) 

hf,u = h^u - JZT^K Vfiu , (2.10) 



yields the expression for the metric perturbation in terms of the newtonian potential (|2.6|) 



00 



d-3^N _ d-3 sr 1 



"o» = -F^i^ = ''rE 



„=_oo (r2 + (z + nL)2) 2 

hji = ; 6in = — /loo ^ii > (2.11) 

^ d-2 iTT ^ d-3 ^ ^ ^ 

where the Latin indices stand for the spatial components, 
^or more precisely po to the power gcd{2,d — 3). 



10 



In the 5d case one can express the newtonian potential as 



$ 



G^M 



sinh (^) 



Lr cosh (^) - cos (^) 



(2.12) 



In figure [^ we give the equipotential surfaces of the newtonian potential in 5d, and they 
look qualitatively the same in any dimension d > 5. 

We close this section with two com- 
ments. First, at higher orders in the per- 
turbation procedure the form of the equa- 
tions is dominated by the linearized equa- 
tions and is given by 



Og 



(m),fiu 



F{9,dg). 




»-r 



Figure 3: Equipotential surfaces of the newtonian 
potential in 5d. 



where [m) is the order under study 
and F(g,dg) are source terms which are 
quadratic, at least, in lower order met- 
ric components and their derivatives. The 
second remark is that in this section all 
the expressions for the metric components were in the harmonic gauge. In the next sec- 
tions we use the Schwarzschild gauge in the near zone. To avoid cluttering the notation 
we will not introduce always a different letter for every gauge — in subsection iA. we give 
different notation for p in the two gauges (p in the Schwarzschild coordinates and ph in the 
harmonic gauge) but later we omit the difference in the notation. However, it is important 
to remember the difference in the gauge between the two zones. 



3. Black hole perturbations 

As explained in the introduction the zeroth order in the near (horizon) zone is the d- 
dimensional Schwarzschild-Tangherlini metric |25] 



ds' 



f{p)dt^ + j^^dp' + p^d^l_^, 



where /(p) = 1 



d-3 

Po 



^d~2 



37:3, Po is related to M via (|2.7D and 



d^ti^o = dx^ + sin^ x dOf -|- • • • -|- (sin^ x sin^ 6*1 .. . sin^ ^d-4) d6d^3 ' 



is the metric on S" . 

In this section we find the linear static perturbations for the d-dimensional 
Schwarzschild solution. Regge and Wheeler pOl derived the linear equations that describe 
small perturbations to the four dimensional Schwarzschild black hole, and here we general- 
ize their method for the static case. These perturbations can be interpreted as deviations 
of the black hole from spherical symmetry due to remote masses. In the case of a compact 
dimension, we are interested in the influence of the black hole (images) on itself. Therefore 
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we assume that the SO(d — 2) symmetry in the spherical coordinates (^i, . . . ,^^-3) is stih 
preserved and the deformation of the black hole takes place only in the (p, x) plane (which 
in "cylindrical" coordinates will be part of the (r, z) plane). We denote in this section the 
d-dimensional Schwarzschild metric by gfj,^ and the perturbation metric by /i^j^. Therefore, 
hnu is a function only of p and x- 



The linearized vacuum Einstein equations can be brought to the simplified form [40| 



SR 



111/ 



\D^D,h - ^D'^D^h^, + D'^Df^.h^y = , 



where h 



g^'^h^^, D^ is the covariant derivative with respect to the background metric 



9fii 



and A 



{al3) 



{Aaf} + Apa) /2 stands for the symmetric part of a tensor A^p- 



d-2 



3.1 Spherical harmonics on S 

Our goal is to simplify the equations by reducing them to a system of ordinary differential 
equations. Following Regge and Wheeler we start by expanding the solution, h^j^u-, into 
generalized spherical harmonics on the sphere S'^~^. Each component of h^i, is transformed 
under local coordinate changes of S'^^'^ like a scalar, a vector or a tensor. Hence, we 
decompose /i^jy into 3 types of spherical harmonics: scalar, vector and tensor harmonics. 
The scalar components are: hu, hpp and hpt 



{hp^,h^ 



pei 



.hn 



The vectors are: 
The tensor is formed of the block 



{htx,htei,--- ,hted-3) ^^^ 



(K 



* 



h^0^ 

^0102 
^6*2 6*2 



h, 



\ 



\ * * * ■■■ f^ed-sdd-s I 
where we denote by * symmetric components. Counting degrees of freedom we find that 
we should have 3 elements in the basis of the scalar harmonics, 2 (d — 2) in the basis of the 
vector harmonics and 2 ™ '^'^ basis of the tensor harmonics. Together we obtain 

all the — ^ — - components of hpy. 

Since we are interested in expansion to harmonics which are static and have SO(d — 2) 

are: 



symmetry the non- vanishing components of the tensor hp,y 



\ 



hpp hpx 



h 



XX 



\ 



he. 



(3.2) 



"d-'i'^d-'i / 

where in addition, the SO(d — 2) symmetry implies that the only independent angular 
components on the diagonal are h^x and hg-^^g^. The rest of the angular components are 
obtained through the relations 

,2, 



ha 



sm 



i-l ha 



,d-3. 
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Note that there are only 5 independent components of h^^, namely hu, hpp, hp-^, h^^, 
ho-^o-^, comprising of 2 scalars, 1 vector and 2 tensors, and these numbers are independent 
of the dimension. Accordingly, we will find 5 linear ordinary differential equations (the 
equations of motion) . 

3.1.1 Scalar harmonics 

By considering the flat {d — 1) -dimensional Laplace equation we can get both the scalar 
spherical harmonics on S*^"^, and the leading radial profile of the multipoles in the nearly 
flat region p ^ po- Since we assume SO{d — 2) symmetry we consider a Laplace equation 
for a function which depends only on one angular variable x- Thus, the Laplace equation 
for a function ^(p, x) on M would be 

dp (/-2 sin('^-3)(x) dp^ip, x)) + p"-' d^ (sin('^-3)(x) d^^{p, x)) = ■ (3.3) 



Separation of variables ^^(p, x) = R{p)^{x) gives us two separated equations for each 
eigenvalue / which denotes also the number of the multipole in the expansion 

oo 

*(p,x) = E^'(^)no''(x). 

1=0 

Hq (x) is the angular function which is associated with each one of the multipoles. 
The extra zero index in the angular function stands to remind us that had we not had the 
SO((i — 2) symmetry, the angular functions, which depends on X) would have additional 
indices (additional "quantum numbers"). (For an example see the 5d case later in this 
subsection). The equation for Hq (x) is 

fl^ + {d-3) cot(x)^ + / (/ + d _ 3) n[,''^ = , (3.4) 

dx ax 

where I {I + d — 3) are the eigenvalues. This equation can be brought to the form of 

a Legendre equation [^] in d — 1 dimensions using the substitution t = cos(x) (then 

the solutions of the equation are called Legendre polynomials of the variable t in d — 1 

dimensions). The solutions can be expressed by a Rodriguez formula (see pl| ]) 



<(X) = ^,^,1 :\, sin^'Hx) -TT f ^^^''^'-\X) , (3.5) 



^(i-^) sin^--^(x)f^-^V 
2'r(/ + f-l) ^""^ Vsm(x) dxj 

where the prefactor with the Gamma functions fixes the usual normalization of the Legendre 
polynomials. 

The functions Ri{p) give us the radial part of the expansion. Ri{p) is obtained as the 
solution of the eigenvalue equation 

±(^p'i-^Ri'{p)'^-l{l + d-3)p''~'Ri{p)=0. (3.6) 

Therefore 

Ri{p) = Jip'^'+''-'^+Kip', (3.7) 

where Ji and Ki are constants. For fixed / and d the first term is the multipole of a mass 
distribution at /? = and the second term is the multipole of a mass distribution at infinity. 
Note that / = is the first multipole — the monopole, / = 1 is the dipole and so on. 
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Since the linear equations for h^j^y ( ^.1| ) are SO((i — 1) invariant we may separate the 
angular variables, and since the perturbed metric depends only on {p,x)i we may expand 
the scalar components hu and hpp in spherical harmonics as follows 

oo 

1=0 

oo 

hpp = J2 Blip) ^\i\x), (3.8) 

1=0 

where the radial functions Ai(p), Bi{p) satisfy some differential equations to be discussed 
later. 

5d scalar spherical harmonics. As a concrete simple example, where we can give 
explicit formulae for the scalar harmonics, let us consider the 5d case. Gerlach and Sengupta 



used this type of decomposition in 5d for the Robertson- Walker spacetime |42|. Let us 
denote the scalar spherical harmonics in 5 dimensions^ by Qf^{x,0,if). They can be 
separated into a product of two types of functions 

QUx,e,^) = u';'\x)Yimi9,v), 

where Yim{0,^) are the usual spherical harmonics on S"^ and II"' (x) ^^^ the "Fock" har- 
monics (see [^, 2^), which are given by 

U^'^ix) = sin'(y) ^'+Mcos((n+l)x)) 
' ^^' ^^' d(cos(x))'+i ■ 

Since we require SO (3) symmetry for the (6, (p) two-sphere, we take Z = m = in the 
spherical harmonics. Thus we arrive to ( |3.4D in the 5d case (with the index n instead of I) 

^!i|! + 2 cot(x) ^ + n (n + 2) uf = , (3.9) 

dx^ dx 

whose solutions are Chebyshev polynomials of the second kind 

rrn,5/ N sin[(n+l)x] _ j. , . ^^ 

^0 U) = r^^ = C/n(cos(x)) . 

sm(x) 

3.1.2 Vector and tensor harmonics 

Given a family of scalar harmonics one can form a family of "scalar derived" vector har- 
monics simply by taking its gradient. There are other, more involved, vector harmonics as 
well, but we shall see now that the "scalar derived" family suffices for our purposes. For a 
concrete example of a basis of vector harmonics in 5d (on S^) see appendix ^ 
Due to the symmetries the vector has a single component ( p.2| ) 



/ip/. = /ipx(P'X) V' (3-10) 



^Note that in 5 dimensions we use different notation for the indices; we keep the index I for the usual 
spherical harmonics Yim{6,tp). Thus, we wiU use fig' instead of Hq' in the general case. 
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where // runs over the angular coordinates only. Denoting h{p,x) = f^ dx' hp^{p,x') we 
have 

hp^ = df,h{p,x)- (3.11) 

Therefore, after expanding h into spherical harmonics and substituting in ( 3. 11]) we 
can write the expansion of hpy. into radial functions Q(/o) as 



hp^ = ^Q{p)dX'ix). 



(3.12) 



A similar argument holds for the tensor components which are essentially /i^^, 
hg-^^Q-^ ([3.2|). Again there are exactly two "scalar derived" tensor harmonics 



t 1 ' 



(3.13) 



and 



nl.d 
.^. ' 



(3.14) 



where ^ij is the metric on the sphere S'^~'^, Di is the covariant derivative on S'^~'^, and the 
second term is proportional to the trace of the first. So we have two tensor harmonics to 
decompose into, which is exactly the number we need, and indeed one can verify that these 
two families always suffice. 

Hence we can decompose the tensor part into radial functions Di{p), Ei(p) 



/ 92 



E^K/') 



/=0 



sin(x)cos(x)^ 



sin(x)cos(x)sin2(6ii)^ 



dx 



\ 



■•■/ 



n[,''(x) + 



/i 



+E^'(/') 



/=0 



sin^(x) 



sin^(x)sin^(6'i) 



V 



\ 



■7 



n[,''(x) 



(3.15) 



The discussion above led us to the conclusion that the only elements of both the vector 
and tensor basis which survived under our symmetry requirements are the "scalar derived" 
ones. These are in the "scalar type" representation under the SO{d — 1) rotation group, 
where by "scalar type" we mean representations whose Dynkin indices are [/,0, ... ,0] for 
some I > 0, namely those in the Z-times traceless symmetric product of the vector repre- 
sentation. More generally, Kodama and Sasaki |^J] pointed out that one can classify the 
different basis elements into three groups with respect to their different representations 
under the rotation isometry group: scalar, vector and tensor "type".^ This classification 



^The vector and tensor types are defined in analogy with the scalar type: a vector type representation 
has Dynkin indices [/, 1, 0, . . . , 0] , namely the traceless product of the 2nd rank antisymmetric representation 
with the i-times traceless symmetric product of the vector representation, and a tensor type is [/, 2, 0, . . . , 0]. 
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into representations of the isometry group sliould not be confused with the classification 
that we used above with respect to local coordinate transformations of 5 . So, according 
to this classification, ( 3.12 ) is a vector of "scalar type" under the rotation group and ( 3.15| ) 



is a tensor of "scalar type" as well. Since we deal with linear equations two different rep- 
resentations cannot mix. Thus even if there were any representations of non-scalar "type" 
in the decomposition of the perturbed metric, they would not appear in the equations for 
the scalar type radial functions Ai, Bi, Ci, Di, Ei. 

3.2 The choice of gauge 

We have now 5 radial fields {Ai{p),..,Ei{p)) defined in ( |3.8[ ), ( p. 121) , ( |3.15| ) for each mode of 
the expansion. We can reduce further the number of fields to 3 using the gauge freedom. 
The gauge transformations of linearized general relativity about a solution g^,y are of the 
form 

where ^^ is an arbitrary vector field — the generator of an infinitesimal transformation. 
The most general generator consistent with the symmetries is 

6 = 0, 

oo 

1=0 

oo 

1=0 

eei =0,...,ee,_3=0. (3.16) 

It is natural to eliminate the "scalar derived" functions Ci{p) and Di(p) in each mode, 
putting h^i, in a diagonal form. We term this the "no derivative gauge". Thus, 

~ "-px + ^P.X + ^X.P ~ "^ ^ PX ^x 



f2(ci + Li + Ml--^Mi^ Ufix) 



1=0 
"■XX ^ "-xx + 2 i^xsx 

— "-xx + ^ ^X,X ~ ^ '^XX ^P 



Y, {Ei + 2pfLi) U'^'^ix) + [Di + 2 Ml) n|)'*'(x) • 



1=0 

The required gauge implies that the functions Li{p) and Mi{p) should satisfy the following 
conditions 

Ci + Li + Ml--Mi = 0, 
P 

Di + 2Mi = 0. (3.17) 
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Moreover, since Ei is arbitrary function of p we can redefine it using the transformation 
El —> El + 2 p f Li. These equations for Li (p) and Mi (p) , although including differential 
operators, are actually algebraic and have a single solution (for each mode) without any 
additional gauge freedom, i.e., the gauge is completely fixed. Applying the gauge transfor- 
mation to any component of (/ig^g^,. . . ■:ho^_.^o^_,^) yields the same equations as in the case 
of h^^ using (|3^). 



3.3 The field equations 

3.3.1 The master equation 

After gauge fixing the perturbed metric diagonalizes and we are left with 3 metric functions 
htt-: hpp, h-^x-' which are functions of {p,x) 

V dx^" dx" = A{p, x) dt^ + B{p, x) dp^ + E{p, x) dQ/_2 



J2 (Mp) dt^ + Blip) dp" + Ei{p)dnl^ li^ix) ■ (3.18) 



It turns out that Einstein's equations simplify if we use h^u instead of h^jy. Thus we 
define A = —-r, B = f B, E = ^ where as usual f{p) = 1 §773. The ansatz in the new 



variables reads 



fMp,x) 



h^„ 



\ 



B{p,x) 



Eip,x) 



(Mp) 



/=0 



E{p,x)---J 



Blip) 



Eiip) 



\ 



jl,d 



K^ix) . 



(3.19) 



Ei{p)...J 
We substitute this ansatz into the linearized Einstein equations, ( |3.lD and using (^ 



we obtain, as expected, five equations 
SRtt = 0; 
2l{l + d-3)Ai-i3p^f + 2{d-2)pf)A'i + p^fB'i-id-2)p^f'E'i-2p^fA'{ = 0, (3.20) 



SRpp = 0; 



21 [l + d - 3) Bi - 3 p^ f A'l + 2 {d - 2) pf B[ + p^ f B[- 



6Rpx = 0; 



A{d-2)pfE'i-{d-2)p'f'E[ + 2p'fA'l-2{d-2)p'fE'{ = ^, (3.21) 
{pf'-2f)Ai + 2{d-3)fBi+pf'Bi-2pfA'i-2{d-3)pfEl = 0, (3.22) 
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6R 



XX 



0; 



wzs ~ ^^xx — 0; 



6Kp)n[,''^(x)-^aKp)nr'(x) = o, 



jl,dl/ / 



sin^(x) 



- i aiip) (cot(x) Ufix) - U'f'ix)) = , 
where we define 

a{p) = A + B + {d-4:)E = tr(/i) - 2E , 

oo 
1=0 

k{p) ^ (d - 3) (i?, - 5,) + 1 p/ (^; - Bi) + \p^r {El + i?n + 

+ {d-2)pfE'i-h{l + d-3)Ei. 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



Three of the equations are second order in the derivatives. These are the evolution 
equations. One of the equations is first order {6Rp-^ = 0) — the constraint equation. For 
/ > the expression in the brackets in ( 3.24| ) is non-zero and we get the following algebraic 
relation 



ai{p) = Ai + Bi + id-4)Ei=0. 



(3.27) 



The case / = is degenerate and we will discuss it separately. 

Now the variables can be separated in ( |3.23| ) and it becomes a second order ordinary 
differential equation 



blip) = (d - 3) (El -Bi) + ^pf{A[- B[) + ^p'f' {E[ + E[') + 
+{d-2) p f E[-U{1 + d-2>)Ei = Q. 



(3.28) 



Using the algebraic relation to eliminate Bi and its derivative from the other equations, 
equations ( 3.20 ) and ( 3.21| ) become 



-ci[p)+l{l + d-3)Ai-p'f'[d-2,)E[ = Q, 
-Clip) -l{l + d-3)Ai-l{l + d-3){d-A)Ei- 
-E'l [{d-2)' pf + p' f {d-3)] - {d-2) p' f E'l' = 0., 



(3.29) 



where we use the abbreviation 

Clip) ^ {2p^ f + {d-2) pf) A'l + p' f A'l' . 
^From the last two equations we can express Ai (p) in terms of Ei (p) and its derivatives 

(3.30) 



A.^J;izM.H^-2)El^,E!]^^^E. 



2l{l + d-3) 
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Now using the remaining equations ( (|3.22| )-( |3.23| )), we arrive to a second order linear 
differential equation for Ei where we introduce x = j-, a dimensionless variable 

ai{x) Ei{x) + pi{x) E'lix) + ^i{x) E'{ {x) = , (3.31) 



where 

ai{x) = -I {I + d - 3) x[{l + d - 2) {I - 1) x'^-'^ - {d - 4) {d - 2)] , 

I3i{x) = (d - 2) (/ + d - 2) (/ - 1) x^-i -f + {d-3)l-2{d- 2f] x^-{d- 2f x^-^ , 

7z(x) = {l-l){l + d-2)x'^+[{d-2)-{l-l){l + d- 2)] x^ -{d- 2) x^~'^ . 

Summarizing, we reduced the equations of static (/ > 1, scalar type) perturbations into 



a single second order ordinary equation for the function Ei (3.31)). This is our "master 
equation" in the variable x. The full solution for the perturbation metric h^u can be 
constructed from the solution for Ei using ( 3.30 ) and (|3.27 ). 



3.3.2 The monopole perturbations 

The solution for / = is equivalent according to Birkhoff's theorem (to leading order) 
with a Schwarzschild black hole possibly with different parameters (mass and asymptotic 
potential) but without a change of its shape (higher multipoles). 

In this case we have IIq' (x) = 1. Hence, the algebraic relation is absent. In this case 
we have only three equations 



6Rtt = 0; 



6R 



pp 



{3p^f' + 2{d-2)pf)A'o + p^f'B'o-{d-2)p^f'E'o-2p^fA'^ = 0., 



3p'f'A'o + 2{d-2)pfB'o + p'f'B'o-A(d-2)pfE'o 



(3.32) 



{d-2)p^f'E'o-2p'fA'^-2{d-2)p^fE'^ 



(3.33) 



6R 



XX 



0; (or 



0) 



{d-3){Eo-Bo) + ^pf{A'o-B'o) + ^p''f'E'o + ^p^fE'^ + id-2)pfE'o = 0. (3.34) 



Subtracting ( |3.32| ) from ( p. 33 ) we obtain a simple equation, where we use again the 
dimensionless variable x = p/po, for later convenience 

A'q + B'o-2E'o-xE'^ = 0. 

The solution of this equation is the following relation between the perturbation metric 
components 



Ao + Bn-En-xE'n 



-2Ci, 



where Ci is a constant. One can verify by substitution that this is the solution of ( 3.32| ) 



as well, i.e., (3.32) does not give us any additional restriction on the solution. Altogether 
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we have 2 independent equations for the three fields Aq, Bq, Eq. Namely, the solution is 
determined up to an arbitrary function. This freedom is exactly a residual gauge freedom 
expected of a Schwarzschild solution p -^ p(/o): it is a result of the unbroken SO(d — 1) 
spherical symmetry. Thus, the general solution should contain, in addition to the free 
function, two constants corresponding to a change in the mass, and an asymptotic potential 
(or equivalently, a rescaling of t). One of the constants will be set by matching, and the 
other one corresponds to a change of parametrization of the branch of solutions, namely, a 
change of "scheme". 

In analogy with the case / > 1, we choose to fix the gauge by the requirement that the 
algebraic condition will be satisfied for the monopole as well, i.e., 

ao{p) = Ao + Bo + {d-4:)Eo = 0. (3.35) 

This eliminates all the gauge freedom apart from a constant shift of the radial coordi- 
nate, which will be a third constant in the general solution. Combining the previous two 
equations together, we obtain a simple differential equation for Eq 

{d - 3) Eq + X E'q - 2Ci = . 

The solution contains two constants which we denote by Ci and C2 

^o{x) = -, — - + 



(i - 3 x'^-s ■ 

Since we use the same algebraic condition as in the I > 1 case, we can obtain the same 
result using ( 3.31| ) for / = 0. Substituting this result into ( |3.343^ ) (and solving it), we 



obtain all the metric components for the monopole perturbation 



,{o),rf _ o r', _ C2 {d - 3) C3 

2 nj.2d—& rj.a 



htt' - 2 Ci - ^ 2d-6 + 7d^ ' (3.37) 






2C.(^i^ + ^V^f^ + ^-^^' '^ 



where C3 is an additional constant. We fix the gauge constant by requiring that the location 
of the horizon remains po, namely that at the horizon the coefficient of the singular second 
order pole (ex 4^) in hpp vanishes. Therefore, we set 

C72 = -^(2Ci + C3). 
a — 3 

Now, in addition, we choose C3 to be 

C3 = -2Ci. 

These choices are advantageous because now one can check that the solution to the 
linearized equations (for the monopole) becomes an exact solution of the Einstein equations 

.<r = ^, (3.39) 
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/iS?'' = 2Ci/, (3.40) 

h^^y = 7^\l ■ (3.41) 



That is, the same form of solution is valid to any order of expansion of Einstein equations 
in the near zone. In each order of the perturbation only one constant, Ci, is needed to be 
determined for the monopole correction and it will be done through the matching process 



to be described in section 4.2. 



3.3.3 The solution of the master equation 

For / > 1 let us make the following change of variable in ( p. 31 ) 

X = x'^'^ , (3.42) 

to obtain a Fuchsian equation with four regular singular points^ 



d^Ei , ^2 , 1 ^__\dEl x + {d-4)wi,d ^ ^ 



where 



and 



wi. 



(3.43) 
d-2 



(/ - 1) (/ + d - 2) 
I 



Pi,d 



d-3 

This is the simplified form of the master equation, from which all metric components can 
be obtained. 

The four regular singular points of the equation are: 0,1, wi^d and cxd. This form of 
Fuchsian equation is known as Heun's equation (see appendix ^ for a review of Fuch- 
sian, Hypergeometric and Heun's equations including the Riemann P-Symbol). It can be 
characterized by the Riemann P-Symbol 

/ 

X, qi, 






1 


Wl4 


oo 











-Pi,d 


-1 





2 


1 + Pl,d 



\ 

where qi^d = -Pi,d {Pi,d + '^) {d - 4) wi^d- 

The exponent difference is an integer number at all the points, except for infinity. Thus, 
the solutions can be represented at each singular point (except for infinity) in two forms of 
power series expansions: one which is just a regular Taylor series (the exponent is zero), 
and a second which may have a logarithmic divergence at the singular point. At infinity, 
the difference of the exponents is an integer when / is a multiple of -^. Then we have 



In the case I = 1 (3.31) becomes an equation with 3 regular singular points (0,1 and oo), namely a 



Hypergeometric equation. The solution of this equation describes a translation of the black hole along the 
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one solution without a log singularity at infinity and second solution which may have a log 
singularity at infinity. Otherwise, the two exponents correspond to two different solutions 
without log at infinity. For physical perturbations of the black hole we will consider only 
the regular solutions at the horizon which correspond to the regular solution at the singular 
point X = 1. 

The characteristic exponents in Heun's equation are determined by 4 parameters of 
the equation (see appendix 0) which we denote by a, (3,^,S . In our case 

a = -Pi,d , 
(3 = 1 +P14, 
7 = 2, 
5 = 1. 

Our case is a rather special case of Heun's equation because its solutions have a very tight 
relation to the solutions of an Hypergeometric equation with the same a, (5, 7 parame- 
ters (the Hypergeometric equation is determined by 3 parameters, see appendix pi). The 
Riemann P-Symbol for this Hypergeometric equation is the following 

1 00 

Fi{x) = p\ -pi^d ; X 

-1 1 1+PM 
The Hypergeometric equation for Fi{X) is 

X(l-X)^^!§^ + 2(l-X)^+p,,(l+p,,) = 0. (3.44) 

Let us define the linear operator C\^^^^^ as 

^Hucn = -HFo + \^ + ^ 7 " ^ 7^ 3^ " Pl,d (1 + Pl,d) 



dX"^ \X X-l X-widJdX '''"^ '^''''' X{X-1){X -wid)' 

(3.45) 

so that ( ^.43 ) can be written as 



'^Heu<n -E'i — , 

and the linear operator CL 

£[lyper ^^(1-^)^+2(1-^)^+ Pl4 (1 + Pl,d) 



SO that ( f3.44 ) can be written as 

^Hyper-^/(^) = 0- 

By defining two additional first order linear operators 



Si^l-id-3)X^, (3.46) 

, ^ ^-3 _^ _ wi^d [d-A) + X 

2-X-ldX X{X-l){X-wi^d)' ^ ' 
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one can verify by direct calculation that 



-^Heun '^1 — ^2 -^Hyper ■ 



(3.48) 



Therefore, if Fi{X) is a solution of the Hyper geometric equation (3.44), then SiFi{X) is a 
solution of Heun's equation ( |3.43| ). Namely, the solutions of ( p. 43 ) can be written as 



Ei{X)=Fi{X)-{d-?,)X 



dFijX) 
dX 



(3.49) 



Notice that the singular point at wi^d does not exist in the Hypergeometric function. It is 
added to the solution through the operator 82- 

The solution regular at the horizon X = 1 written in terms of a hypergeometric 
function is 



E, 



l-{d-3)X 



_d_ 
dX 



(1 - X) 2F1 { 1 
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d-3 
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I 
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.2;1-X 



(3.50) 



and is our final expression for the solution of the master equation ( 3.31 ) . Note that when 
^33 is a natural number the solution for Ei{X) is a polynomial and thus regular at all 3 



singular points. We add here the expression for Ai ( |3.30| ) written in terms of the variable X: 



{d-3)Hd-2){X-l) f dEi d^Ei 

' 2l{l + d-3) \ dX dX^ 



d-4 ^ 



(3.51) 



4. The matching procedure 

4.1 Transformation of coordinates: Sch^varzschild to harmonic 

The choice of harmonic coordinates in the asymptotic zone as opposed to Schwarzschild 
coordinates in the near zone leads us to consider carefully the transformation of coordi- 
nates in order to compare the multipole moments in the overlap region correctly. The 



harmonic condition (2^) is satisfied in the asymptotic flat region only in cartesian coordi- 
nates. Thus, we have to work in coordinates which are cartesian in spatial infinity. Since 
the Schwarzschild coordinates are spherical coordinates at infinity we take the d — 1 carte- 
sian coordinates that are related to (p, x, ^1, • • • , Qd-z) where p is the Schwarzschild radial 
coordinate. 

Now one can verify that by the transformation of the form p — > ph{p) we are able to 
transform the Schwarzschild coordinates to harmonic coordinates, provided 



dp 



d-2 



PPt' 



dph 
dp 



(d-2)/"^p;,. 



(4.1) 



( |4.1| ) is obtained by substituting any of the expressions for x'^ expressed in spherical coor- 
dinates, where we take ph{p) instead of p, into the definition of the harmonic gauge (static 
case) 

Ax^ = 0. (4.2) 
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Without loss of generality we take 

Az = Aiphip)cos{x)) = 0. 

Since the laplacian of any of the cartesian coordinates reduce to ( [4.1D , we know that the 
transformation p — > Ph{p), when ph{p) is defined as a solution of (0), gives the required 
transformation. Equation ([4.1D has two solutions which we can write as a power series. 



When p^ pq, (4.1) becomes 



Substituting ph = p we get a quadratic equation for the characteristic exponents at infinity 
(see appendix (^ ) 

k'^ + k{d-3)-id-2) = 0, 

whose solutions are 

k = l,_(d-2). 

We set one boundary condition that for large p, ph ~ p, i.e., for large p we expect that the 
two gauges coincide. The leading terms in the transformation are 

4.2 Zones and basic dialogue arrows 

A black hole in a spacetime with compact dimension is distorted due to tidal fields that 
originate from its effective mirror images. Namely, the periodic boundary conditions in 
the compact dimension change the structure of the solution compared to the case of fiat 
asymptotic infinity. We assume that the black hole is small compared to the period of the 
compact dimension L (po ^ L). Hence in the limit p ^ pQ -^ L the exact solution should 
become a d-dimensional Schwarzschild solution. For pQ -^ p ^ L^ far from the black hole, 
we have a flat d-dimensional Minkowsky spacetime with a periodic coordinate (which we 
denote by z). 

In the intermediate region we can approximate the exact solution by matching two 
perturbative expansions from both sides. This procedure can be viewed as a dialogue of 
multipoles where the black hole changes its shape (mass multipoles) in response to the field 
(multipoles) created by its periodic "mirrors", and that in turn changes its field and so 
on. Usually when we use an iterative method for solving non-linear equations, we have to 
supply enough data about the boundary conditions at each step of the iteration so that the 
solution would be determined completely ,i.e., without any free parameters or functions. 
Here the two known solutions (the Schwarzschild solution and the Minkowsky spacetime 
with a periodic coordinate) will play the role of boundary conditions for our problem. 
Hence, each perturbative expansion will complete the missing information in the other 
expansion, order by order. After describing the method we will show that it is well-posed, 
namely that all the information needed for a certain perturbation order is provided by the 
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Figure 4: The division of the spacetime into two overlapping zones: the near zone p ^ po where the 
perturbative parameter is L^^, the asymptotic zone p^ <^ p ^ L where the perturbative parameter 
is po and the overlap region in between. 



previous ones. The approximate solution in the end will be written as a combination of 
a post-newtonian expansion far away from the black hole and an expansion of the same 
order around the asymptotically flat black hole solution. For this purpose, we first divide 
the spacetime into two zones with an overlap region (see figure ^ : 

The near zone: for p ^ po we expand the solution around the d-dimensional Schwarz- 
schild solution characterized by po (see section ^. Since p ^ pQ in this zone, the small 
dimensionful parameter of the expansion in this zone would be L~^. Notice that when 
L — > oo (i.e., the limit of "unfolding" the compact dimension) we are left only with the 
first term (zeroth order) of the expansion which is the d-dimensional Schwarzschild solution. 
The next terms represent the influence of the compact dimension on the black hole's metric. 
From the point of view of the near zone, as the period L is large compared to po, the terms 
in the expansion series become smaller and the convergence of the perturbation series is 
better. 

The asymptotic zone: for po <C p ~ -L we expand around the flat background with 
a compact dimension of size L. Here we take pQ the size of the black hole as the small 
dimensionful parameter (this is the so called post-newtonian expansion, see section 0, with 
the matching added as a boundary condition). The terms in the expansion represent the 
distortion of the flat spacetime (with a compact dimension) due to the presence of the 
black hole and its "mirror images" . From the point of view of the asymptotic zone, as the 
mass of the black hole becomes smaller, the terms in the expansion series become smaller 
and the convergence of the perturbation series is better. 

The overlap region: each one of the expansions cannot be determined independently 
since it requires boundary conditions at the overlap. Hence, it carries free parameters that 
should be determined by matching to the other expansion in the overlap region po <^ p <^ L. 
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This is achieved by comparing the multipoles of the expansions order by order; In the 
asymptotic zone we look at the multipoles in the limit p <^ L and match them with the 
multipoles of the near zone expansion in the limit p^ pQ. 

There are three parameters that determine the matching procedure for each term in 
the multipole expansion: 

• The order of the multipole expansion which we denote by n. A term of order n in 
the asymptotic zone expansion is proportional to Pq while a term of the same order 
in the near zone is proportional to (l/L)". 

• The multipole number /. In each order a given multipole can have contributions from 
lower multipoles (via non-linear sources from the same zone, e.g. addition of angular 
momenta) and from matching which is linear and as such can be performed only 
between identical multipoles in both zones. 

• The number of dimensions d. 

The interplay between these three parameters l,n,d creates the matching pattern which 
we describe below. 

In the previous sections we obtained the linear perturbations around the exact solutions 
(equations (|3.43|) and (2.11)). In order to get non-linear corrections by the method of 



successive approximations we substitute in Einstein's equations expansions of the form 

where g^iJ is the exact solution. For each term of this expansion we get an equation for 
the metric components of the form 

A5&)) =-^n(5("-^),5("-'\...,5(^^), (4.5) 

where the superscript n denotes the order in the perturbation, C is linearized operator 
(which is second order and independent of n) and J-'n stands for nonlinear terms that are 
functions of the previous orders and their derivatives. Since the asymptotic zone expansion 
is in the harmonic gauge, C there is simply the laplacian. In the near zone, as we saw in 
section |3|, the operator is more involved and it is given by >C^g^jj, defined in ( |3.45D . At each 
order one is free to add to the metric a solution to the homogeneous equation (an element 
of ker(£)). The homogeneous solution of any order n can be expanded into multipoles. 
Thus we can write the general form of the n-th correction to the metric in the form 

oo 

# = Ci,) + E -^^ f!^^'(p) no''(x) . (4.6) 

1=0 

The first term is the particular solution of order n. The second term is the multipole 
expansion of the homogeneous solution where IIq (x) is the spherical harmonics function, 
defined in ( |3.4D (remember that we have SO{d — 2) symmetry). /^" (p) is the multipole 
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(n) 

radial function and c^ are coefficients that we should determine using the matching pro- 
cedure for each multipole I and order n (in d dimensions). For the matching we would have 
to expand the particular solutions to multipoles as well. 

iFioui the previous subsection we learn that in the limit p ^ po the Schwarzschild 
and harmonic coordinates coincide. Hence we are allowed to compare leading terms in 
the two zones without coordinate transformations. Moreover, in this limit (p ^ po) the 
homogeneous linear operator in the near zone becomes the laplacian (Einstein's equations 
in this limit coincides with the newtonian poisson equations). Therefore, in both zones, 
the leading terms will be multipole functions of the laplacian, namely, leading terms that 
their dependence on p goes like p' or ^"('+"^"3) ^ggg (^)). 

In the homogeneous solutions of the near zone the leading terms (the ones to be 
matched) behave like p^ (when one takes the limit p ^ po). The component with the 
subleading behavior p~''+"~3) is determined by the boundary condition of regularity at the 



horizon. This can be seen explicitly from the solution that we have for Ei in section 3.3.3. 
The Riemann P-Symbol (appendix ^) shows us that the asymptotic behavior of the two 
power series expansions'^ at infinity is />' or p^C+'^^^j^ p^,, ^j^g asymptotic zone the role is 
reversed - the leading term in the overlap zone {p <C L) is p^('+'^^'^) while the subleading 
term p' is fixed by the boundary conditions at infinity. 

We summarize the discussion above using a graphic representation of the multipole 
expansion for multipole number I in figure ^. The numbers stand for the order in the 
multipole expansion. The two rows of cells represent the expansion in the asymptotic zone 
(above) and in the near zone (below). The left side of each cell contains the leading term 
in the overlap region (that should "receive" the information from the matching) . The right 
side of each cell contains the appropriate term to match the left side of different cell. The 
horizontal part in the cell reminds us the order of perturbation — pg in the asymptotic 
zone and L"" in the near zone. The matching procedure actually becomes a process of 
determining which cell's left side receives information from which cell's right side (as we 
will see, not all the cells are relevant for this process). 

To determine the pattern of perturbation, namely the basic "dialogue arrows" of in- 
formation flow in a diagram such as figure |5|, the basic observation comes from dimensional 
analysis. A mass multipole moment of multipole number /, Aii, is defined as 



Ml, .-/-..... 



where the indices ij denote the cartesian coordinates, and the superscript is a reminder 
that these are multipoles at the origin. The gravitational field of multipole of order / of 
mass distribution at the origin behaves like 



pl+d-3 ' 



and the multipole moment itself has the length dimensions 

[A^O] =/ + (i-3. (4.7) 



^Here for simplicity of the matching procedure we consider the equation 3.43 with dimensions, namely, 
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Figure 5: The representation of tiie expansion in the near and asymptotic zones for a multipole 
moment I, before the matching. Each cell in the upper row stands for a term in the asymptotic 
zone expansion and each cell in the lower row stands for a term in the near zone expansion. Each 
cell consists of the small parameter in the appropriate power and two terms to be matched. 

Such mass multipoles at the origin are used to supply matching boundary conditions for 
the asymptotic zone. Suppose we are at order n in the asymptotic zone so that 






(4.8) 



where m is the near zone order from which it was matched. Comparing (4/7) and ( [4. 81) we 
conclude that 

n-m= [M'}] =l + d-3, (4.9) 

namely, upward pointing "dialogue arrows" have a step-size oil + d — 3. In a similar manner 
the multipole moments at infinity (used to supply boundary conditions for the near zone) 
yield 

n-m= [Mf°] = -I , 

since the (dimensionless) gravitational field behaves like 

Namely, downward pointing "dialogue arrows" have a step-size of +1. 

Alternatively, the same result can be described from the following point of view. In the 
perturbation process we expand a dimensionless quantity which satisfied Laplace's equation 
in the two zones and compare the expansions in the overlap zone. In the matching process 
we identify two terms that make one dimensionless quantity in the overlap region. The 
process is demonstrated in figure |^. Looking at a cell of order m in the near zone, we 
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Figure 6: Two steps in tlie matcliing procedure. A term of order m and niultipole number I in tlie 
near zone determines the multipole of order n — m + l + d — 3 in the asymptotic zone. Afterwards 
the latter is used to determine a term of order m'~n + l — m + 2l + d— 3in the near zone. 



see that the corresponding term in the expansion is proportional to 



T m ^L-\-d — 3 ' 



In order to 



create a dimensionless quantity, we have to match this term to order n = l + m + d — 3 in 
the asymptotic zone. Then we have a leading term proportional to 

l+m+d—3 

Po 



up to a numerical constant. Using the same argument, this cell in the asymptotic zone 
is matched to the cell m' = n + I = m + 21 + d — 3 in order to create the dimensionless 
quantity 

Pop' 

Apparently, this process can continue to create an infinite chain of matchings for each 
multipole which consists of two unequal steps; from the near to asymptotic zone we have 
a step of An = I + d — 3 orders and from the asymptotic to the near we have a step of 
An = / orders. 

4.3 Pattern of the dialogue 

Let us look in detail into the pattern of the dialogue. We would like to summarize all the 
"information flow arrows" between the various orders and to conclude how the computation 
should be performed to determine the metric up to a prescribed order in one of the zones 
(following what we term "the critical route"). 

First we note that since our problem is defined with reflection symmetry z — > — z, one 
should consider only multipoles that satisfy this symmetry. Applying the transformation 
X — >■ X + TT in (|3.5| ) shows that only even multipoles are relevant for our expansions. 

We summarize the information flow arrows for a diagram such as figure |5|: 
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• Within a zone. 

These are the ordinary non-hnear sources of the perturbation procedure. At order 
n there can be a non-zero source only if there are some non-zero metrics at orders 
Hi < n such that some hnear combination of rii with integer positive coefficients 
yields n. 

• Inter-zone arrows. 

We saw above that arrows going "down" (from the asymptotic to the near zone) 
advance by I orders, while those going in the reverse direction advance by / -|- d — 3 
orders. Since / can only be even there are two basic steps: 2 and d — 3. 

Finally there is an "initial condition" : the leading correction to the zeroth order metrics 
is the newtonian potential at order d — 3 in the asymptotic zone. Therefore the larger d is 
there will be more orders where the metric vanishes. 

^From a practical point of view it is especially important to determine the "critical 
route". Namely, suppose we wish to compute the metric up to a prescribed order n in 
a certain zone. Which orders must be determined on the way? Clearly we should know 
all the lower orders in that zone.^ In addition we should have information from the other 
zone. If the original zone in which we were interested is the near zone we should know 
the metric at the asymptotic zone at orders n — / for all even /, and similarly if it were 
the asymptotic one we need to know all the n — (d — 3) — I orders of the near zone. We 
see that in both cases the limiting I is the monopole I = which sets the critical route to 
arrive from the other zone: to reach the near (asymptotic) order n the critical route passes 
through the asymptotic (near) order n (n — d + 3). It is clear that the procedure does not 
run into closed loops since the sum of the two step-sizes, being d — 3, is greater than 0, and 
in that sense there are sufficient boundary conditions to determine the expansion, namely 
it is well-posed. 

In order to make this analysis more concrete we shall first outline the first steps in 
the procedure for d = 5 which is somewhat special and then do the same for some general 
d>5. 

In the case of 5 dimensions d — 3 = 2 and this implies that only the even orders 
are involved in the expansion (the same holds true for any odd dimension). Hence in 5 
dimensions we actually expand in the small parameters Pq and j-j (see figure 0) . This is a 
great simplification to the perturbation scheme. Moreover, all the even multipole functions 
(see ( 3.50D ) are polynomials in 5 dimensions. 



In figure |^ we show the first steps in the matching procedure in 5 dimensions, until 
the fourth order. In order to have the metric until the fourth order in the two zones we 
need, besides the particular solutions, to match the monopole (/ = 0) and the quadrupole 
(/ = 2). In order to obtain higher orders we will have to add higher multipoles to the 
matching. Note that we have to know the fourth order in the asymptotic zone in order to 
get the fourth order in the near zone. This situation is rather different for d > 5. 



* Actually, this may not be necessary since for example orders n' such that n — (d — 3) < n' < n will not 
contribute. Yet, normally, we are not interested in a correction before all previous ones were determined. 
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Figure 7: The first steps in the matching procedure in 5 dimensions. Until the fourth order only 
monopole and quadrupole matches appear in the two zones. Each one of the matches is represented 
by a different color. The newtonian approximation enters in the second order in the asymptotic 
zone. It is clear from the figure that the monopole is the limiting and hence the most "fresh" 
multipole that should be determined for a certain order. 

In 5d it is easy to understand that the method is well-posed. If we look at a cell of order 
n in the near zone, then all the (even) niultipoles from I = to I = n are determined from 
the previous orders in the asymptotic zone when the monopole is the last to be matched 
— it is determined from the same order but in the asymptotic zone. We do not have to 
determine any multipole higher than n because they do not appear at this order — in order 
to create a dimensionless quantity out of p"+'^ for A; > 1 we have to multiply it at least by 
,„_).fc (thus this multipole will appear at order n + A; at least in the near zone). Therefore, 
for any order until n the metric in the near zone is determined. Similar argument is valid 
for the asymptotic zone as well, let us look at a cell of order n in the asymptotic zone. 
All the multipoles from / = to / = n — 2 are determined by matching to the near zone 
when / = (the monopole) is the last one to be matched (it is matched with the cell of 2 
orders before) — see figure |^. Therefore, the information for any cell comes from the cells 
of lower orders in the other zone in 5d. 

The first steps in the matching for d > 5 dimensions are shown in figure ^ Since the 
newtonian approximation enters at the order d—3 > 2, the post-newtonian correction comes 
at 2 (d — 3) . The first quadrupole contribution in the near zone comes in (d — 3)-|-2 = d—1. 
Then, as opposed to the case of d = 5, the second matching of the monopole in the 
near zone comes after the quadrupole matching, since 2{d — 3) > d — 1. As we go to 
higher dimensions, there are more multipoles in the near zone that we can match with 
the multipoles in the newtonian approximation before adding the second correction to the 



31 



Newtonian 
Approximation 



Monopolc Match (t=0) 
Quadrupolc Match (1=2) 



d-3 




Asymptotic 2(d-3) 
Zone 




d-3 



d-2 



d-l 



Near 
Zone 



Figure 8: The first steps in the matching procedure for d > 5. Compare with the previous figure 
(the same in 5 dimensions) . Here the first quadrupole matching in the near zone is done before the 
second correction to the monopole. 

monopole. This can be calculated precisely: The length of a step of a match from the 
asymptotic zone to the near zone is An = /. The distance between the newtonian and 
post-newtonian approximations is An = d — 3. So any multipole that satisfy I < d — 3 
comes before the second correction to the monopole. An interesting limit is when d — > oo. 
Then, the newtonian approximation is enough to determine all the orders of expansion in 
the near zone. 

5. Matching results 

After describing the procedure in general, let us start and calculate the first coefficients of 
the approximate solution. Note that in the matching process it is enough to match a single 
metric component to guarantee the matching of the whole metric. For our purposes it will 
be convenient to consider for the matching the dimensionless quantity g^^./p'^. 

5.1 The first monopole match 

Actually there is a zeroth monopole match which relates the mass measured asymptotically 
(through the newtonian potential) with the po parameter of the Schwarzschild solution. 
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namely a matching between the zeroth order in the near zone and order n = d — 3 in the 
asymptotic zone (the first red arrow in figures ^ and ^) . This match was done in section ^ 
(based on ||) and gives the identification p^'^ = ^3^0^- 

The first monopole match is between the monopole in the newtonian approximation 
(n = d — 3) and the first monopole correction to the d-dimensional Schwarzschild solution. 
Being the / = correction the spherical symmetry of the black hole will not be lost, but 
rather the physical interpretation is that the black hole is placed in a region with a constant 
non-zero gravitational potential due to the array of image black holes. 

First we expand the newtonian approximation^ in multipoles around the origin includ- 
ing not only the monopole term but also the quadruple in anticipation of the computation 
of eccentricity. The expansion around p = in polar coordinates (p, x) is 
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where C, is Riemann's zeta function^'' and 
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(5.1) 
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according to the Rodrigues formula (^I^). The superscript "(asymp), d — 3" indicates that 
this is the correction of order d — 3 in the asymptotic zone. In the same way we will denote 
the expansions in the near zone by the superscript "(near)". Note that there is no dipole 
term in the expansion (or any odd multipole) due to the symmetry X ^ X + '^• 

Recall the expressions that we had for the metric monopole perturbations ( 3.39| )"( |3.41| ). 
Comparison with the monopole term in ( |5.1| ) gives us the value of the matching constant 
Ci for the first correction to the metric in the near zone (n = d — 3) 
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The final form of the metric in the near zone with the first monopole correction is 
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(5.2) 



^Recall that p in the asymptotic zone is in harmonic gauge while p in the near zone is in the "no 
derivative" gauge around the Schwarzschild gauge. Since we compare only leading terms in the two zones 
and in order to simplify the notation we denote both of them by p. 
^"Riemann's zeta function is defined as C{s) = X]^i ^- 
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As mentioned above this monopole correction can be attributed to the non-zero new- 
tonian potential at infinity of the near zone. The same result can be obtained in a less 
formal way by rescaling of t reflecting in this way the change in the asymptotic boundary 
conditions. 

Since there is no change in the shape of the black hole we can measure the effect of 
the correction by calculating the correction to the dimensionless quantity k Ad-2 where 
K is the surface gravity of the horizon and Ad-2 is its {d — 2)-dimensional area. 

To find the area we note that the horizon in this approximation is located at a constant 
Ph which is defined as the root of guip)=0. Using the approximated metric above we get 
that the horizon does not move (in our gauge) 



P^ 



Pt' + O 



1 
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The area of the horizon is 



a—z 
Ad-2 = 9xxiPH) 2 ^d~2 

-PO ^^-2^1 + (^ _ 3) ^.-3 + ^ 

The surface gravity calculated using (|B.1D of appendix O is 
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Substituting the approximate metric we get 



P=PH 



d-3 / {d-2)C{d-5) pt' 
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Ld-i 



Combining the results together we obtain 




(5.5) 



5.2 Matching multipoles to the newtonian approximation 

Our next goal is to obtain an expression for the deviation from spherical symmetry — the 
eccentricity of the horizon. For this purpose we have to match the quadrupole term in 
the near zone to the newtonian approximation. In this subsection we give a more general 
recipe, matching any multipole to the newtonian approximation. 

As was explained in the previous section, this matching is important for higher multi- 
poles as we increase the number of dimensions. From the matching scheme we learn that 
the multipoles in the newtonian approximation match with An = / orders higher in the 
near zone expansion. 

Let us count the number of constants to match. Given / there are two radial solutions 
in the near zone. Regularity at the horizon fixes a particular combination and so it remains 
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to set the normalization of this combination by comparing the leading terms from both 
zones. Thus we need to determine a single constant per multipole. 

In the near zone we have the multipole functions written in terms of solutions of a 
Hypergeometric equation (recall ( p. 44 ) and ( 3.49| )). Regularity of the metric at the horizon 



requires that the relevant solution of the Hypergeometric equation be regular there. The 
horizon is located at X = 1 when we use the variable X as we defined in (3.42) 



Using the Riemann P-Symbol of ( 3.44| ) we can obtain the regular solution at X = 1 
(see Q) and substitute it into ( p. 49] ) 

i.,= (l-(d-3)xA)(i_x)^,(l--^,2 + ^,2;l-x), (5,6) 

The correction to the near zone metric of order I + d — 3 is proportional to Ei 

(near),/+d-3 
^-2^^ = A pi Ei{X) Ili\x) , (5.7) 

where Di is the matching constant which depends on the multipole and the dimension. 
This constant will be determined by matching to the asymptotic zone. If d — 3 > / the 
newtonian approximation suffices and otherwise non-linear terms from lower orders should 
be added. 

In order to match this solution in the overlap zone we have to find its leading term 
as X ^ oo. The regular solution of the Hypergeometric equation (13.44 ) at X = 1 can be 



written as a combination of the two solutions at infinity, using the identity (see |^ 
2Fi{a,l3,r,l-X) = ci{X-iy"2Fi (a,a-j + l,a-p + l;j^^ + 



+C2 (X - 1)/Fi T/?,/? - 7 + 1,/? - a + 1; ^ 



X 

where 

r(7)r(/3-a 



(5.8) 



ci 



C2 



r(7-a)r(/?)' 
r(7)r(a-/3) 



r(7-/?)r(a)- 

Hence, in our case the leading term of the Hypergeometric function is 

r(2)r(/\ + i) , , 

r(i + ^)r(2 + ^) 

Substituting into (|5.6| ) we obtain the leading term in Ei as X ^ oo 

r(:^ + i) ^^^' 



Ei- ^.. . T'U , . (^-1)(-) • (5.9) 



r(i + ^)r(2 + ^)^ Vpo. 

Note that the exponent ~ p of the leading term is exactly what we anticipate for multipole 
function of order I. The determination of the prefactor allows one to obtain the constants 
Di by matching with the asymptotic zone. 
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I 




Figure 9: The effective "mirror" images induce tidal forces that stretch tlic horizon of the black 
hole. We quantify this by the eccentricity which is defined using two sections of the horizon. 

5.3 The eccentricity 

The lowest order deviation from spherical symmetry of the near zone metric occurs in 
n = (d — 3) + 2 = d— 1, where we match the quadrupole (see figures ^ and ^). We can 
quantify deviation of the black-hole from spherical symmetry by the eccentricity of the 
horizon. We define the eccentricity of the horizon as in [^ 

e = ^-l, (5.10) 

where Au is the (d — 3)-dimensional area of the equatorial sphere at x = f and A± is the 
{d — 3)-dimensional area of the polar sphere at x = Oi ^ (see figure P) . Since g^^^^ is a 
function only of p and x (recall the symmetry SO(d — 3)) we obtain 

^11 = ^d-3 9xx [P = PH,X = ^ 

A^ = n^-A / 5xx (P = PH)sm ~^(x)dx- 
Jo 
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Matching the quadrupole in the near zone in order n = d — 1 to the newtonian approxi- 
mation n = d — 3 (equating the quadrupole terms — 1 = 2 — in { f).7\ ) and ( |5.1| ) using the 
leading term ( |5.9D of the quadrupole radial function in the near zone) we obtain 



D, 



(d-2)c(d-i)r(i + ^)r(2 + ^) p^ 



d-l 



n^s + i) 



L'^- 



Therefore 



(near),d-i 



9xx 



2A, 



p- 



D2E2{X)Ii'^\x)- 



(5.11) 



(5.12) 



For the calculation of the eccentricity we need the value of E2{X) at the horizon, namely, 
X = 1. From ( |5.6D we can see that 



E2{X = l) = d-3. 



(5.13) 



As the first non-spherical term is at order n = d — 1 in the near zone, the leading term 
in the eccentricity is proportional to ,'^_i . Thus, we have to substitute g^^ i^ the 

d-l 

expression for the eccentricity ( p. 101 ) and find the coefficient of j^^rj in the expansion. First 
we expand 



d-3 

9xx 



P=Ph 



24, 



pV [l + --- + ptD2{d-3)Iil'\x) + 



d-3 
1 \\ — 

1? 



Pr (i + --- + ^V^pg^2nf(x) + o(^ 



and at X = "^/^ we get 



d-3 

9xx 



P=PH,X=2 



"^'i'---T0/«"^-°{h 



(5.14) 



(5.15) 



where the ellipsis stand for monopole corrections which, of course, do not contribute to the 
eccentricity. Substituting in ( 5.10| ) we obtain 



id-3f 



pID2 



n 



d-4 



2 ''^ ^ V^d-3 Jo 

The integral gives us 



d-4/ 



ni'{x)sin'-'{x)dx + 



,_„+0[§,]^ (5.16) 



T2,d/ 



U',^^{x)sin''-''{x)dx = - 
d 



1 (d-3\ (2^^Hd-l)n^ dV^ \ 



This using identities of the Gamma function (see P^ j) yields the final formula for the 
eccentricity 

(5.18) 



(d- 


- 3)^ (r(2 


+ 3^))'C(d- 


-1) 


Po 
Ld-1 


+ 0(#)- 




8{d- 


^)n^,) 
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Figure 10: A plot of the function ei(d) which gives the d-dependence of the first coefficient of 
proportionahty in the expansion of the eccentricity, norniahzed as in (5.20). Note that ei(d) is 
monotonously decreasing to 0. 



We can also write the eccentricity as a function of a differently normalized variable: 



the relative size of the horizon in conformal coordinates which we define as (see |20|) 



2ph 2d-3po 



L 



L 



(5.19) 



where ph is the location of the horizon in conformal coordinates ds^ = e^ dt^ + e^ ■^ (dr^ + 
which were used extensively in the literature and will be used here only 



dz^) + e 



2C 



dnl. 



for setting the normalization of rj. rj is the same quantity denoted by x in |20, 21|. One 
gets 



(d-3)Mr(2 + ^))2c(d-i) 



(d-4) (d+1) 
d-3 



(^-2)r(;^) 



-^ + 0(7]'' 



eiT] 



^d-1 



(5.20) 



The constants ei{d) are illustrated in figure IC and in the next table 



d 


ei 


5 
6 


fi ^ 2.886 
~ 2.099 


7 
8 


2100 - 1-^^^ 

~ 0.947 


9 


~ 0.603 


10 


~ 0.375 
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5.4 The "Archimedes" effect 

We define the inter-polar distance -Lpoies to be the proper distance between the "poles" of 
the black hole measured around the compact dimension 

r-L/2 

-lpoies = 2 / dz^/g^;, (5.21) 

Jzh 

where zh denotes the location of the horizon. It is convenient to define a dimensionless 
parameter out of Lpoies 

y ^ 1 _ b^ . (5.22) 

Thus y is the relative decrease in the size of the compact dimension at r = due to the 
presence of the black hole. Like any other property of the system it is a function of the 
single dimensionless quantity which characterizes it, which we choose to be ij, the relative 
size of the horizon in conformal coordinates defined in ( |5.19| ) . 

y{rj) is an interesting quantity to compute in our framework since it involves both 
zones in an essential way. The way to compute it is to pick some mid-point Z, divide the 
integration between the two zones 

I Vies = r dz a/^F^ + r dz A^F^ , (5.23) 

and confirm that the result is independent of the choice of mid-point (up to the specified 
order in r]). 

Clearly y{j] = 0) = 0. Here we wish to compute y to first order in r/. Which order 
is required on each zone for the evaluation? Near the horizon it is enough to take the 
zeroth order, that is the Schwarzschild solution with no corrections (a factor of r\ multiplies 
it as all distances scale with /jq and the definition of y includes a division by V). In the 
asymptotic zone one would need in principle the first order correction (in r\)^ but since the 
leading order newtonian potential is of order r/ and d > 5 it suffices to consider again 
the zeroth order, namely flat compactified space. 

Having the matching with the asymptotic zone in mind we write the Schwarzschild 
metric in the following conformal coordinates 

i_ ^^2 (near) ^ _ {^_^ ^^2 + (^ + ^)^ (^^2 ^ ^2 ^j^2_^) (524) 



where 



P 
These coordinates approach the newtonian gauge as p ^ ph- 



^=[Pli]'\ (5.25) 



We may now compute the contribution from the near zone to (5.23) 



1 /" , . / (near) 



, , dz\/gTr '= (5.26) 

^ -'zH 
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transforming to a dimensionless (" = z/ p^ 



1 \ <*-'* 



C' 



d-3 



dC 



(5.27) 



adding (+1 — 1) to the integrand in order to separate the part which diverges with Z 

(5.^ 



Ph 
L 



Z \ fZ/ph 

1) + 

Ph 



c 



d-3 



1 + T^ 1 '"' - 1 1 dC 



transforming to tf; = 1/C in order to facihtate taking the hmit ph ^ rj ^ 



^,Ph 

L L 



Phiz 



l + y^'^-sV-^.i 



dw 



■w^ 



Z Ph J 



where the definite (and finite) integral Id is 




(5.29) 



(5.30) 



and at this order we neglect the remainder of the integral 



Ph 
' L 



Ph/Z 






The contribution of the (flat metric) asymptotic zone to ( |5.23D is simply 

LJz ^^'^ "2 L' 



(5.31) 



(5.32) 



where due to the low order of the calculation and the choice of coordinates we can use 
the same value for Z in both patches, with no need for corrections arising from a further 
matching of the patches. 

Summing ( 5.29 ), ( 5.32|) according to (^.23|) and using the definitions ( ^.22]) , (|5.19| ) we 
confirm that the mid-point dependence, (Z/L), drops (up to the relevant order) and we 
get our result 

y = 2IdV + o{v), (5.33) 

where the function Id is shown in figure |l^. It is monotonously increasing starting from 
I5 = 0, Iq = .6845, tending to 1 as d — > 00. The special behavior at 5d where the "amount 
of space" outside the black hole remains fixed to first order in the black hole size was 



termed "the black hole Archimedes effect" [20[ since the black hole seems to "repel" as 
much space as its own size (figure [ll| ) . 

In principle we could have used the corrections to the metric computed in section g and 
subsection |5.2| to improve on this computation and determine the next to leading correction 
as well. 
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Figure 11: Illustration of the "Archimedes effect". The small black hole repels an amount of space 
proportional to its size. In 5 dimensions, the black hole repels (to first order) an amount of space 
which is equal to its size (measured in conformal coordinates). 



0.8 



0.4 



0.2 




Figure 12: The function Id (5.30,5.33) representing the decrease in the inter-polar distance (to 
first order) relative to the size of the small black hole (in conformal coordinates). Note that /a = 
and it rises monotonously to 1. 
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A. Hypergeometric and Heun's Equations 



Let us consider a linear ordinary differential equation of second order 

(P f df 



(A.l) 
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where a{x), (3{x),'y{x) are holomorphic functions. When a(xo) = and /3(xo) 7^ (or 
7(xo) 7^ 0) then x = xq is a singular point of ( [A.l| ). This point is called regular singular 
point if the limits 

,. (X - Xq) P{x) 

hm = Ai , 

x^xQ a{x) 

and 

(x-xo)^7(x) 
hm = A2 , 

x^xo a{x) 

exist. The regular singular point at infinity is defined in the same way for t = in the 
equation 

"2 1 /3(i) 



— - + 



t t2 „(1) 



dt^t^a{\y 



which is obtained after changing the variable x into t = -. Equation ( [A.l| ) is called a 
Fuchsian equation if all the singular points of the equation are regular. 
The roots of the equation, 

^i'^ + {Kl-l)^l + K2 = Q, 

are called the characteristic exponents at xq. There are two linearly independent solutions 
to the equation. Thus, at any singular regular point, there are two different expansions 
into power series. If the difference of the exponents is not an integer, the power series near 
a singular point xq are of the form 

00 



(X - Xo)'' ^ fflj (x - Xo) 



i=0 

where /i is one of the characteristic exponents. When the difference of the exponents is an 
integer one of the expansions may involve a logarithmic term. 

Any Fuchsian equation with 3 singular points can be transformed, by transformations 
of the dependent and independent variables, to a standard form where the singularities are 
at 0, 1 and 00 

x(l-x)— 4 + [7-(a + /3 + l)a;]3^-a/3/ = 0, (A.2) 

ax^ ax 

where a, /3, 7 are parameters. The equation in this form is called the Hypergeometric 
equation. The solutions of ( [A.l| ) can be characterized completely by the singularities and 
the corresponding characteristic exponents at each singularity. There exists a compact 
scheme to summarize the information about the solutions of the Hypergeometric equation 
— the Riemann P-Symbol (see for instance ||45| , ^) 



P Q ; X 








1 




00 












Q 


-7 


7- 


-a- 


-P 


P 
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The first row of the matrix indicates the three regular singular points. The two numbers 
beneath each singular point are the characteristic exponents at each singular point. The 
expansion around of the solution with a zero exponent is the Hyper geometric function 

where (a)„ = Y^(a) '^hen a is different from a negative integer. The other 5 expansions 
can be obtained by transformations of this Hyper geometric function. 

A Fuchsian equation with 4 singular points can be transformed similarly to a canonical 
form where the singularities are at 0, 1, a and 00 (a can be any point in the complex plane 
different from the other singular points). In this canonical form the equation is called 



Heun's equation ^, ^ 

(P f (1 ^ l + a + /3 — 7 — (5\ df aBx — q „ ,. , 

ax^ \x x — l x — a J ax x [x — 1) [x — a) 

This type of equation is a generalization of the Hyper geometric equation. Then in a similar 
manner a Riemann P-Symbol can be constructed for the 4 singular points: 

/ 1 a cx) \ 

P a ; X, q . 

yi-7 l-S j + 6-a-P (3 y 

Unlike the Hyper geometric equation, Heun's equation is not characterized completely 
by its characteristic exponents at each singular point. There is another parameter — q. 
This parameter is called the accessory or auxiliary parameter. 

B. The surface gravity 

Let us compute the surface gravity for a metric of the form 

ds'^ = -f{p) dt^ + g{p)dp^ + aij dx' dx^ , 



^ = -t{D^^C)iD^Cuj, 



(where aij dx^ dx^ stands for the rest of the spatial part of the metric). Using the definition 

of surface gravity [|^ 

1 

■2' 
where C,^ is the Killing vector field dt- Using the fact that 

r* = ^1— 
^P 2 / ' 

f 

""25' 

we get the following formula for the surface gravity 
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C. 5d scalar harmonics 

For concreteness, let us write down a basis for vector harmonics in 5 dimensions. First, we 
have a natural basis element for vector harmonics which is derived from the scalar spherical 
harmonics 

namely D^ is the covariant derivative on S^ 

This basis element is a vector, whose direction is defined by the gradient on the sphere. 
We look for another two vector basis elements on S^ which are orthogonal to this one.^^ 
We fix the arbitrariness in the choice of the other two vectors by considering the S0(3) 
symmetry in the {6, ip) coordinates. Namely, we choose a second vector basis element which 
is orthogonal to the first one in the {9, ip) plane 

where 

sm V -^ 
V^ = sin(x) e" Uu = (0, ^—r D^, - sin 6 sin x Dq) . 

Here e^p represent the two dimensional Levi-Civita tensor for the coordinates — ip and zero 
for components with x coordinate. The multiplication by sin(x) is due to the difference 
between the inner products on 5^ and S*^. 

The third basis element is defined as a vector which is orthogonal to the other two 
basis elements 

where 

and e^j^up is the three-dimensional Levi-Civita tensor. (In the calculation we used the 
eigenvalue equation for the spherical harmonics on 5^.) Gerlach and Sengupta |^^ obtained 
exactly this form of a basis of vector harmonics on S^ using a different method. 
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